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1.  INTRODUCTION 


Project  “Internal  Tidal  Hydrodynamics  and  Ambient  Characteristics  of  the  Adriatic 
(ITHACA)”  concentrates  on  diurnal  internal  tides  -  their  generation  at  the  Adriatic 
shelf  break,  propagation  across  the  Palagruza  Sill  and  dependence  on  stratification 
and  currents  in  the  area.  Systematic  oceanographic  measurements  started  in  1950s  at 
the  Split-Gargano  transect  by  the  Institute  of  Oceanography  and  Fisheries  (Split)  and 
in  1960s  in  a  wider  area  by  the  Hydrographic  Institute  (Split).  It  was  found  that  the 
area  is  under  the  strong  influence  of  seasonal  and  interannual  variability.  Early  short¬ 
term  current  measurements  revealed  a  pronounced  seasonal  signal,  with  oppositely 
directed  alongshore  flow  in  the  winter  and  summer  seasons  ( Zore-Armanda ,  1968). 
Seasonally  dependent  circulation  contributes  to  a  strong  temporal  variability  of 
thermohaline  properties  ( Grbec  and  Morovic ,  1997).  Since  the  beginning  of  1980s 
long-term  current  measurements  were  occasionally  performed,  showing  barotropic- 
like  current  reversals  on  scales  of  a  few  days  ( Leder  et  al .,  1996).  More  recently, 
higher-frequency  phenomena  attracted  some  interest.  In  particular,  it  has  been  found 
that  during  the  stratified  season,  diurnal  internal  tides  are  generated  by  the  interaction 
of  diurnal  barotropic  tides  with  topography  and  that  they  result  in  thermocline 
variability  with  amplitudes  reaching  3  m  ( Mihanovic ,  2005).  No  indications  were 
found  of  semidiurnal  internal  tides,  probably  because  the  semidiurnal  tidal  currents 
have  a  minimum  at  the  shelf  break.  As  the  detection  of  internal  tides  was  based  solely 
on  data  collected  by  three  thermistors  on  the  island  of  Lastovo,  and  there  were  no 
current  measurements  to  support  the  finding,  an  intention  of  the  present  project  is  to 
provide  a  more  complete  data  set,  which  would  enable  propagation  characteristics  of 
internal  waves  to  be  determined,  not  just  their  amplitudes.  Attention  is  also  paid  to  the 
way  changes  of  background  stratification  and  currents  modify  internal  waves,  and 
vice  versa  -  to  a  possible  influence  of  internal  waves  on  deductions  based  on 
measurements  that  are  scattered  in  space  and  time.  Thus,  observations  and  modeling 
are  planned  so  as  to  enable  the  following  hypotheses  to  be  tested: 

•  Diurnal  surface  tides  generate  internal  tidal  waves  at  the  Adriatic  shelf  break. 

•  Internal  tidal  energy  radiates  away  from  the  shelf  break.  Some  waves  travel 
towards  the  northwest,  thus  crossing  the  Palagruza  Sill. 

•  Generation  and  propagation  of  baroclinic  waves  is  influenced  by  stratification  and 
currents  prevailing  in  the  area. 

The  project  started  on  12  April  2005  with  the  financial  support  of  the  USA  Office  of 
Naval  Research  (award  No.  N00014-05-1-0698)  and  the  Croatian  Ministry  of  Science, 
Education  and  Sports  (contract  No.  3/2005)  and  is  expected  to  last  until  30  September 
2007.  Preparation  of  the  field  work  is  planned  for  the  first  year,  execution  of  the 
experiment  for  the  second  year,  and  analysis  of  the  data  and  model  results  for  the  third 
year. 

This  is  report  on  the  first  year  of  work  on  the  project,  which  -  following  the  project 
proposal  -  concentrated  on  the  preparations  for  the  field  work  by  (l)  analyzing 
historical  data  collected  in  the  Palagruza  Sill  area  and  (2)  developing  a  simple  model 
that  approximates  the  Adriatic  shelf  break  with  a  depth  discontinuity  and  basin 
stratification  with  two  layers. 
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2.  CLIMATOLOGICAL  CONDITIONS  IN  THE  MIDDLE  ADRIATIC 


2.  1.  Background 

The  earliest  oceanographic  research  of  the  Adriatic  Sea  started  already  in  the  19th 
century.  However,  systematic  and  regular  measurements  in  the  middle  Adriatic  began  in 
1950s  by  the  Institute  of  Oceanography  and  Fisheries  in  Split.  One  of  the  most 
investigated  areas  in  the  middle  Adriatic  is  the  Split-Gargano  transect  (Figure  1),  region 
characterized  by  a  strong  temporal  variability  of  thermohaline  structure  determined  by  a 
number  of  factors:  interaction  between  the  atmosphere  and  sea  leading  e.  g.  to  deep 
water  formation  in  the  north  Adriatic,  river  discharge,  water  exchange  through  the 
Otranto  Strait  and  topographic  effect  of  the  Palagruza  Sill  (Zore-Armanda  and  Bone . 
1987). 

Generally,  current  outflows  from  the  Adriatic  to  the  Mediterranean  in  the  bottom  layer, 
whereas  the  Mediterranean  water  enters  the  Adriatic  in  the  intermediate  layer.  In  winter, 
in  the  north  Adriatic,  very  cold  dense  water  is  formed,  whereupon  it  sinks  to  deep  layers 
of  the  Jabuka  Pit  and  is  occasionally  advected  across  the  Palagruza  Sill  (Zore-Armanda, 
1963).  The  transect  area  is  also  under  the  influence  of  saltier  water  advected  from  the 
south  Adriatic  and  from  the  Mediterranean.  Advection  of  the  Mediterranean  water, 
called  "digression"  ( Buljan ,  1974),  brings  Mediterranean  saltier  water  into  the  Adriatic, 
causing  a  strong  increase  of  salinity  in  the  middle  Adriatic.  Since  the  temperature  of 
Levantine  Intermediate  Water  (LIW)  is  higher  than  that  of  the  Adriatic  water, 
"ingressions"  are  also  manifested  as  an  increase  of  temperature  (Zore-Armanda,  1969a). 
Additionally,  ingressions  cause  increase  in  nutrients  and  consequently  enhance  primary 
production.  Relating  these  phenomena  to  the  climate,  Zore-Armanda  (1969b)  stated  that 
the  most  important  factor  enhancing  the  water  exchange  between  the  two  basins  is  the 
horizontal  air  pressure  gradient  over  the  east  Mediterranean.  The  locations  of  the  Iceland 
cyclone  and  Siberian  anticyclone  centers  were  found  to  control  the  air  pressure 
differences  over  the  Adriatic  (Zore-Armanda,  1969c).  It  was  found  that  such  changes 
could  be  related  to  conditions  in  a  wider  area  of  the  north  Atlantic  and  Europe  (Grbec  et 
al.,  2003). 

The  Palagruza  Sill  transect  is  exposed  to  a  number  of  different  impacts,  which  makes  it 
suitable  for  investigation  of  specific  thermohaline  features  on  both  the  seasonal  and 
interannual  scales  (e.g.  Orlic  et  al .,  1992;  Grbec  and  Morovic ,  1997;  Grbec ,  1998;  Zore- 
Armanda  et  aL,  1999). 

In  the  following,  atmospheric  forcing  in  the  middle  Adriatic  is  briefly  considered. 
Hydrographic  conditions  are  then  described  by  presenting  thermohaline  features 
observed  along  the  Palagruza  transect,  which  also  comprises  Stoncica  -  the  station  best 
investigated  in  the  region. 
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Figure  1.  Map  of  the  Adriatic  showing  bathymetry  of  the  basin  and  position  of 
permanent  meteorological  stations  and  stations  CJ07  -  CJ12  of  the  Split-Gargano 
transect. 


2.  2.  Atmospheric  forcing 


The  middle  Adriatic  is  located  in  the  area  in  which  strong  orographic  effects  of 
Dalmatian  mountains  modify  weather  conditions  whereas  the  climate  has  almost 
Mediterranean  features.  Strong  cyclonic  activity  over  the  Adriatic,  especially  in  the  cold 
part  of  year,  causes  different  weather  types  to  occur  (Penzar,  1967;  Loncar  and  Bajic, 
1994)  with  frequent  exchange  of  Sirocco  and  Bora  events  (Figure  2).  During  summer, 
this  area  is  under  the  influence  of  subtropical  high-pressure  belt,  so  warm  and  dry 
conditions  prevail  with  well-developed  sea-land  breeze  systems.  Seasonal  variations  of 
main  atmospheric  parameters  such  as  solar  radiation,  wind,  air  temperature  and 
precipitation  cause  strong  seasonal  variability  of  heat  and  water  flux  between  the 
atmosphere  and  the  sea  (Figures  3-5). 


Figure  2.  Annual  wind  roses  (relative  frequency  of  wind  directions,  left,  and  average 
wind  velocities,  right)  according  to  wind  measurements  at  the  Split-Marjan  station 
over  the  1961-1990  interval. 
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The  annual  cycle  of  heat  flux  is  controlled  by  heat  gain  due  to  incoming  short-wave 
radiation  and  heat  losses  through  outgoing  long-wave  radiation,  evaporation  and 
conduction  of  heat.  These  components  of  heat  flux  for  the  middle  Adriatic  were 
estimated  using  different  semi-empirical  equations  ( Grbec  et  al.,  1997)  using  data 
originating  from  surface  meteorological  stations  distributed  along  the  east  Adriatic 
coast  {Grbec,  1998;  Matic,  2005).  After  several  automatic  meteorological  stations 
were  established,  short-wave  radiation  and  long-wave  backward  radiation  were 
measured  directly  (ADRJCOSM  reports,  2002-2003). 


Figure  3.  Seasonal  cycle  of  upward,  downward  and  net  heat  flux  obtained  from 
measurements  performed  at  Split-Marjan  over  the  1961-1990  interval. 


Month 


Figure  4.  Seasonal  cycles  of  surface  air  and  sea  temperature  at  the  meteorological 
station  Split-Marjan,  for  the  1961-1990  interval.  Vertical  bars  represent  standard 
deviations. 
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Figure  5.  Seasonal  cycle  of  evaporation  for  various  stations  along  the  east  Adriatic  coast 
(lop)  and  precipitation  at  the  meteorological  station  Split-Marjan  (bottom),  for  the  1961- 
1990  interval.  Vertical  bars  represent  standard  deviations. 


Mean  seasonal  net  heat  flux  cycle  (Figure  3)  shows  that  from  April  to  September  the 
sea  receives  heat  from  the  atmosphere,  reaching  maximum  in  July  (134.6  Wm"2  at 
station  Split-Marjan).  The  sea  loses  heat  most  strongly  in  December  (1 12.1  Wm"2  at 
station  Split-Marjan);  at  the  annual  scale  the  loss  of  heat  is  caused  by  long-wave 
radiation  (60%),  evaporation  (33%)  and  much  less  by  conduction  (7%). 
Correspondingly,  from  May  until  August  the  sea  is  colder  and  from  October  until 
March  the  sea  is  warmer  than  the  overlaying  atmosphere  (Figure  4). 
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Water  exchange  at  the  air-sea  interface,  determined  using  monthly  mean  values  from 
meteorological  stations  Dubrovnik,  Hvar  and  Split,  is  shown  in  Figure  5.  Evaporation 
at  different  stations  along  the  coast  is  considerably  different  from  that  at  the  open  sea, 
both  in  annual  courses  and  in  absolute  values.  Evaporation  along  the  east  Adriatic 
coast  reaches  maximum  in  cooling  season,  with  the  values  far  from  the  coast 
(meteorological  station  Hvar)  being  significantly  lower  than  those  close  to  the  coast. 
Fresh  water  input  by  precipitation  decreases  while  going  towards  the  open  sea,  since 
there  is  no  orographic  effect  which  is  considerable  near  the  coast.  Therefore,  due  to 
various  influences,  mean  seasonal  P-E  cycles  differ  significantly  for  the  coastal  and 
the  open  middle  Adriatic  area. 

It  is  interesting  to  note  that  during  extreme  weather  conditions  (Table  1),  for  example 
during  strong  Bura  events,  the  heat  loss  from  the  sea  is  extremely  high,  causing 
cooling  of  the  sea  surface  and  gradual  salinity  increase  due  to  strong  evaporation. 

In  addition  to  the  strong  seasonal  cycles  in  the  atmospheric  and  boundary  layer 
parameters,  long-term  increasing/decreasing  trends  can  be  detected  which  highlight  the 
changeable  middle  Adriatic  climate.  These  changes  are  clearly  reflected  in  the  long-term 
variability  of  hydrographic  properties  as,  for  example,  sea  surface  temperature  (Figure 
6). 


Table  1.  Mean  daily  values  of  heat  flux  components  (W/m2)  during  several  extreme 
weather  events  at  the  meteorological  station  Split-Marjan.  Os  denotes  shortwave 
radiation,  Oh  is  sensible  heat  flux,  Qe  is  latent  heat  flux,  Qi  is  longwave  heat  flux,  and 
Qnet  is  net  heat  flux  at  the  air-sea  interface.  TA  and  Ts,  in  °C,  denote  surface  air-  and  sea- 
temperature,  respectively. 
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Figure  6.  Long-term  sea  surface  temperature  (SST)  fluctuations  obtained  from  twice- 
daily  measurements  at  the  Split-Marjan  station. 


2.  3.  Seasonal  and  interannual  variability  of  temperature  and  salinity 

2.  3.  L  Split-Gargano  transect 

Thermohaline  variations  over  the  year  at  the  Split-Gargano  transect  are  a  result  of 
different  processes  in  the  surface  and  deeper  layers,  causing  presence  of  particular 
water  masses  in  different  seasons,  as  seen  from  T-S  diagrams  (Figure  7).  Different 
therrnohaline  features  observed  from  May  1998  to  May  1999  at  the  Split-Gargano 
transect  are  analyzed  and  compared  with  average  water  mass  characteristics  for  the 
region. 

The  figures  presented  describe  therrnohaline  properties  in  the  middle  Adriatic  in  1998 
and  1999.  In  winter  1999,  enhanced  heat  loss  and  low  Po  River  discharge  resulted  in 
the  North  Adriatic  Dense  Water  (NAdDW)  formation,  in  contrast  to  winter  1998 
when,  due  to  weak  air-sea  exchange  and  moderate  Po  River  discharge,  fresh  water 
was  present  in  the  whole  north  Adriatic  and  consequently  the  NAdDW  was  not 
formed  ( Zore-Armanda  et  al.,  1999;  Vilibic  et  ai,  2004).  After  the  formation, 
NAdDW  flowed  as  a  high  density  current  along  the  west  Adriatic  slope  towards  the 
Jabuka  Pit.  Its  densest  core  interacted  with  the  topography  and  rapidly  sank  when 
reaching  relatively  steep  slope,  west  of  the  Jabuka  Pit.  The  second  part  followed  the 
slope  as  a  high  density  current,  but  not  interacting  with  the  bottom  topography, 
detouring  the  pit  along  its  southeast  side  and  hitting  the  bottom  somewhere  in  the  east 
segment  of  the  pit.  Changes  of  therrnohaline  properties  along  the  Palagruza  Sill  are 
documented  in  Figures  8  and  9,  using  data  collected  between  April  2004  and  February 
2005  during  cruises  of  RV  Bios  (project  Jadran).  Close  to  the  east  Adriatic  coast  less 
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saline  water  is  present  in  the  surface  layer  during  a  greater  part  of  the  year.  This  is 
especially  visible  during  maximum  freshwater  inflow  in  spring.  During  winter 
vertically  homogenous  water  column  is  found  in  the  whole  transect,  with  temperature 
and  salinity  being  lower  close  to  the  east  coast.  Of  course,  themohaline  conditions 
observed  between  April  2004  and  February  2005  can  be  discussed  only  as  a  case 
study,  taking  into  account  prevailing  atmospheric  processes. 
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Figure  7.  T-S  diagrams  constructed  from  data  collected  at  the  Split-Gargano  transect 
between  May  1998  and  May  1999. 
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Figure  8.  Distribution  of  temperature  and  salinity  along  the  Palagruza  Sill  transect  for 
April,  June,  July  and  September  2004.  Gargano  is  on  the  left  and  Split  on  the  right. 
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Figure  9.  As  in  Figure  8,  except  for  December  2004  and  February  2005. 


2.  3.  2.  Stoncica  station 

Long  series  of  monthly  values  of  temperature  and  salinity,  measured  between  January 
1961  and  December  1995,  were  used  to  analyze  mean  seasonal  thermohaline  features  of 
the  Vis  island  area.  The  data  were  collected  during  regular  monthly  cruises  at  the 
oceanographic  station  Stoncica  (see  Figure  1).  Measurements  were  made  mostly  once  a 
month  (only  exceptionally  several  times  per  month)  at  standard  oceanographic  depths 
(0,  10,20,30,50,  75  and  100  m). 

The  mean  seasonal  cycle  of  temperature  (Figure  10)  is  very  pronounced  down  to  20 
meter  depth,  with  rather  small  temperature  changes  at  100  meters  where  temperature 
does  not  fall  below  12°C.  The  curves  show  isothermal  conditions  at  all  depths  in  late 
autumn  and  winter,  whereas  stratification  is  present  from  May  until  October.  Annual 
temperature  maximum  in  the  deeper  layers  is  lagging  the  surface  one,  whereas  minimum 
at  all  depths  occurs  within  the  same  month.  An  earlier  study  ( Grbec ,  1998)  showed  that 
at  the  depth  of  1 0  meters  the  temperature  maximum  has  a  delay  of  9  days  with  respect  to 
that  in  the  surface  layer,  whereas  at  20  meters  the  maximum  is  delayed  48  days.  The  75 
and  100  meter  depths  have  a  lag  of  93  days  relative  to  the  surface.  Minimum 
temperature  in  the  surface  layer  occurs  in  February,  and  at  75  and  100  meters  at  the 
beginning  of  March.  An  interesting  feature,  related  to  temperature  variations  in  a 
particular  month  (represented  here  with  the  standard  deviation),  may  be  observed:  the 
largest  variability  during  summer  and  early  autumn  is  seen  in  the  intermediate  layer. 


12 


(°C)  f  o  T« rr^« Utufn  f*  C)  Tamp«r«l»jr«  (°  C> 


L#r*i  {p- 10)  m 


Lfy*r  (10-20)  m 


Liytt  (20-30)  m 


Uy«r  (30-50)  m 


Month 


Month 


Ljr*f  (50-75)  m 


Uy»r  (75-80)  m 


20 


•  M*jn 


M«jiu30 


24 


20 


*  M«J'i  MuntSD 


24 


22 


22 


L»f*«  (ftS-tOOi  m 


Sljt.an 


a« 


24 


*  Utin 


UunlSD 


22 


20 


10 

14 


*2 


12  3  4  5  0  7  8 


8  10  11  12 


Month 


Manth 


Figure  10.  Mean  seasonal  temperature  cycles  for  selected  layers  at  Stoncica  station. 
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Figure  11.  Mean  seasonal  salinity  cycles  for  selected  layers  at  Stoncica  station. 


14 


Layer (m)  Layer (m) 


Temperature  f  C)-o- 

12  13  14  15  16  17  13  19  20  21  22  23  24  25  26 


Salinity 

38.0  38.2  36.4  36.6  36.8  37.0  37.2  37.4  37.6  37.8  38.0  33.2  38.4  38.6  38.8  39.0 


0-10 

10-20 

20-30 

30-50 

50-75 

75-85 

85-100 


Temperature  C)-o- 

12  13  14  15  16  17  18  19  20  21  22  23  24  25  26 


S  a  li  n  ity  — • — 

36.0  38.2  36.4  36.6  36.3  37.0  37.2  37.4  37.6  37.8  33.0  38.2  38.4  38.6  38.8  39.0 


0-10 

10-20 

20-30 

30-50 

50-75 

75-35 

85-100 


Figure  12.  Mean  vertical  temperature  and  salinity  profiles  for  Stoncica  station. 
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Figure  13.  Time  evolution  of  temperature  and  salinity  obtained  from  monthly  mean 
values  for  the  Stoncica  station  and  the  1961-1995  interval.  Depth  layers  -1  to  -7  denote 
layers  having  0-10,  10-20,  20-30,  30-50,  50-75,  75-85  and  85-100  meter  depths, 
respectively. 
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The  seasonal  salinity  cycle  (Figure  11)  shows  salinity  increase  with  depth.  Again,  the 
largest  seasonal  variations  are  observed  in  the  surface  layer,  whereas  at  greater  depths 
salinity  is  more  steady.  In  the  surface  layer  down  to  10  meters  there  are  two  minima,  in 
May  and  in  July.  Salinity  maximum  is  not  pronounced,  although  larger  values  at  all 
depths  are  observed  in  late  summer  and  early  autumn.  The  observed  interannual 
variability  is  probably  related  to  advection  from  the  Mediterranean. 

The  vertical  distribution  of  temperature  shows  the  presence  of  thermocline  below  the  20 
m  depth,  while  the  distribution  of  salinity  shows  halocline  at  almost  the  same  depth 
(Figures  12  and  13). 

In  the  deeper  layers  the  maximum  temperature  values  in  May  are  the  consequence  of 
fast  transport  of  heat  to  the  layers.  In  summer,  thermocline  suppresses  transport  of  heat 
to  the  deeper  layers  and,  consequently,  lower  temperatures  are  observed  below  the 
thermocline  throughout  the  summer  (Figure  12).  Besides  the  heat  flux,  the  wind  has 
considerable  influence  on  temperature  of  the  water  column.  Increase  of  wind  speed  in 
autumn  coincides  with  the  temperature  increase  in  deeper  layers,  as  wind  induced 
mixing  transports  heat  from  the  warm  surface  layers  to  the  deeper  ones. 

Salinity  in  the  surface  layer  during  the  heating  season  mirrors  P-E  values,  whereas  - 
somewhat  surprisingly  -  in  autumn  and  winter  salinity  and  P-E  values  are  in  phase. 
When  the  thermocline  is  well  developed,  lighter  east  Adriatic  surface  waters  reside  in 
the  surface  layer  and  are  advected  to  the  Split-Gargano  transect.  They  bring  about  a 
decrease  of  salinity  in  the  surface  layer,  mixing  being  prevented  by  the  fully  developed 
thermocline.  In  summer,  due  to  current  system  in  the  intermediate  layer  ( Zore-Armanda , 
1969b),  advection  of  saltier  water  from  the  south  Adriatic  (Mediterranean)  is  observed 
(Figure  12).  In  the  warmer  part  of  the  year  the  two  salinity  minima  may  be  attributed  to 
the  Neretva  River  outflow  (in  May)  and  to  the  influence  of  the  rivers  outflowing  from 
the  south  Adriatic  coast  (in  July). 

It  is  evident  that  vertical  mixing  processes  have  essential  influence  on  thermohaline 
properties  of  the  whole  water  column,  but  with  different  strength  over  the  year.  These 
vertical  processes  are  to  a  considerable  extent  controlled  by  direct  atmospheric  forcing. 
In  addition,  thermohaline  properties  are  heavily  influenced  by  horizontal  advection, 
which  is  determined  by  oceanographic  as  well  as  meteorological  conditions. 
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3.  ANALYTICAL  MODELLING  OF  INTERNAL  TIDES 
3.1.  Two-layer  analytical  model  with  a  step-like  topography 


The  first  simple  analytical  model  describing  the  generation  of  internal  tides  was 
presented  by  Rattray  (1960).  He  described  the  analysis  of  two-layer  system  with  an 
abrupt  change  of  topography  under  presumption  of  long  waves,  and  showed  that  this 
model  gives  an  amplification  of  internal  tides  with  respect  to  surface  tides. 
Introduction  of  real  topography  and  continuous  density  distribution  limits  the 
possibilities  of  analytical  modeling,  making  numerical  approach  the  only  possible 
alternative.  Analytical  solutions  of  idealized  model  were  also  considered  for  linear 
shelf  slope  between  the  coast  and  depth  discontinuity  (Rattray,  1960),  and  for  linear 
discontinuity  slope  ( Prinsenberg  and  Rattray ,  1975).  Modal  approach  described  by 
Prinsenberg  and  Rattray  (1975)  is  appropriate  for  linear  shelf  slope  and  any  given 
stratification,  whilst  Baines  (1982)  presented  analytical  models  that  are  applying 
different  approach  in  internal  tide  research,  based  on  the  ray  theory.  Standard  integral 
equation  is  solved,  and  solutions  can  be  obtained  for  any  topography  but  only  for 
certain  (continuous)  density  distributions.  The  greater  part  of  models,  analytical  and 
numerical,  was  vertically  two-dimensional  (Maze,  1987;  Sherwin  and  Taylor ,  1990), 
but  Baines  (1986)  and  Sherwin  et  ai  (2002)  emphasize  importance  of  3D  modeling, 
because  of  the  influence  of  along-slope  barotropic  flow.  Obliqueness  of  barotropic 
wave  was  considered  by  Guizien  et  ai  (1999).  They  were  following  original  Rattray’s 
approach,  with  different  boundary  conditions  at  the  coast.  Namely,  Rattray  used  wave 
reflection  at  the  coast,  which  gave  partial  standing  waves  on  the  shelf,  while  Guizien 
et  al.  in  their  analysis  of  internal  tide  generation  by  surface  tidal  wave  incoming  to  the 
shelf  with  certain  incidence  angle,  used  radiation  boundary  condition  at  the  coast  (no 
reflection  at  the  solid  boundary). 

Analytical  approach  to  the  generation  of  internal  tides  in  the  vicinity  of  the  shelf  edge 
is  relevant  if  certain  assumptions  are  fulfilled.  These  assumptions  are  two-layer  fluid 
and  rapid  depth  change  (Guizien  et  al. ,  1999).  Long  waves  and  step-like  topography 
assumptions  are  valid  as  long  as  wave  lengths  are  significantly  larger  than  shelf 
dimensions  along  rapid  depth  changes.  The  island  of  Lastovo,  from  which  internal 
tides  in  the  Adriatic  were  observed  for  the  first  time  (Mihanovic,  2005),  is  located 
about  140  km  from  the  center  of  the  South  Adriatic  Pit,  with  an  abrupt  depth  change 
at  a  distance  of  about  30  km.  Since  wave  lengths  of  barotropic  surface  waves  are  of 
the  order  of  1000  km,  it  can  be  concluded  that  this  assumption  is  satisfied.  Another 
important  supposition,  which  is  very  important  for  this  kind  of  analytic  models,  is 
two-layer  fluid  approximation.  It  also  implies  the  existence  of  the  first  mode  of 
internal  oscillations  only,  which  is  the  case  for  the  best  part  of  shallow  seas,  where 
stratification  exists  when  seasonal  thermocline  is  developed  ( Huthnance ,  1989). 
During  June  2001  the  most  significant  temperature  oscillations  were  detected  at  the 
depth  of  22  m  (especially  during  the  first  half  of  June),  and  oscillations  at  15  and  36 
meters  were  not  that  significant  (Mihanovic,  2005).  It  could  be  concluded  that  in  the 
interval  mentioned  thermocline  was  positioned  in  this  part  of  the  water  column.  CTD 
measurements  from  August  1990  confirm  typical  vertical  distribution  of  thermohaline 
properties  of  this  area,  with  three  typical  layers.  However,  thermocline  (pycnocline) 
was  very  pronounced  and  narrow,  so  vertical  structure  can  be  approximated  by  two- 
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layer  system  ( Leder ,  2002).  Therefore,  this  assumption  in  a  wider  Lastovo  area  is  also 
fulfilled. 


Two-layer  analytical  model  used  in  theoretical  study  of  the  generation  of  internal  tides 
in  the  Adriatic  is  based  on  decomposed  equations  for  barotropic  and  baroclinic  mode 
of  oscillations  of  a  two-layer  fluid  in  rotating  rectangular  channel  (e.g.  Csanady, 
1971;  Orlic,  1988).  Traditional  and  hydrostatic  approximations  were  introduced, 
whilst  tidal  forcing,  nonlinear  effects,  lateral  friction  and  stress  at  the  pycnocline  were 
neglected: 
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where  £ p ,  up  and  v p  represent  elevation  and  velocity  components  for  the  surface 
(barotropic)  mode,  u'u  and  v'  elevation  and  velocity  components  for  the  internal 
(baroclinic)  mode,  p  and  h  are  constant  density  and  undisturbed  depth  of  the  surface 
layer  (thermocline  depth),  p  and  h'  are  constant  density  and  undisturbed  depth  of 
the  bottom  layer,  rxs.  and  t  are  surface  stress  components,  and  t'xb  and  t'yb  are 
bottom  stress  components.  Notation  e  =  {pl—  p)j p  is  introduced  as  well. 

The  connection  between  solution  for  the  surface  mode  and  surface  elevation  and 
current  component  oscillations  in  the  bottom  layer  is  given  by: 


Ur=ltr’  Vr=V?' 
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while  relation  between  the  internal  mode  solutions  and  surface  elevation  and  current 
component  oscillations  in  the  surface  layer  is  given  by: 


r  n  h'  ,  ti  , 

=  °;  «»  =  —ruu >  vu  -  rv« 
h  h 


Integral  solution  for  the  surface  and  bottom  layer  can  be  reconstructed  by  combining 
solutions  of  (3.1)  and  (3.2): 


u=up  +Uu 
V  =  vp+vu 
u'  =  u'p+u'u 

v  =vp  • 


(3.3) 


Having  in  mind  that  theoretical  study  of  internal  tides  in  the  Adriatic  will  be  based  on 
frictionless  two-layer  analytical  model,  variations  in  direction  parallel  to  the  shelf 
discontinuity'  can  be  neglected,  together  with  the  surface  and  bottom  friction.  Coriolis 
force  components  were  usually  disregarded  in  traditional  studies  of  conversion  of 
barotropic  into  baroclinic  signal  (Miles,  1974),  which  is  valid  for  the  surface  mode 
since  Rossby  radius  of  deformation  is  one  order  of  magnitude  larger  than  spatial 
dimensions  in  which  generation  processes  take  place  (Gill,  1982).  Specifically, 
internal  tides  are  generated  at  the  edge  of  deeper  Adriatic  region,  in  which  surface 
Rossby  radius  of  deformation  has  an  order  of  magnitude  of  about  1000  km,  while 
generation  processes  take  place  in  area  of  about  100  x  100  km.  This  argument  is  not 
valid  for  internal  mode,  since  in  June  2001  internal  Rossby  radius  of  deformation  near 
Lastovo  island  was  about  3  km  ( Mihanovic ,  2005),  and  lengths  of  internal  tidal  waves 
are  usually  between  20  and  40  km  (e.g.  Holloway,  1983,  1987;  New  and  Da  Silva, 
2002).  For  this  reason  Coriolis  force  has  an  important  influence  on  prevention  of 
soliton  generation  (solitons  are  nonlinear  waves  connected  to  internal  tides,  Gerkema, 
1996).  However,  since  analytical  modeling  described  in  this  report  is  aimed  at 
understanding  the  basics  of  internal  tides  in  the  Adriatic,  Coriolis  force  will  be 
discarded  for  both  the  surface  and  internal  modes.  Equations  (3.1)  and  (3.2)  then 
become: 
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surface  mode: 


(3.4) 


internal  mode: 
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Elevation  and  velocity  in  the  surface  and  bottom  layer  are: 

h  +  h 
ti  , 

U  =  Up+UrJ  =Up-  —  Uu 

ti  =  up+u'u  =  UP+K- 


General  expression  for  surface  mode  elevation  will  be  £ p  —  Ap  exp|/(&tf  ±  kpx)\,  and 
for  internal  mode  elevation  =  A'uexp[i(cot±k'ux)\.  By  combining  equations  (3.4) 
and  (3.5),  wave  equation  for  the  surface  and  internal  mode  is  obtained: 


surface  mode: 


internal  mode: 


g[h  +  ti) 
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a2<r: ..  hti  a2c 
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which  for  the  assumed  wave  forms  give  connection  between  wave  number  and 
angular  frequency: 


(3.5) 
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surface  mode: 


internal  mode: 


P  g{h  +  h') 

,2  _  {h  +  //) 

“  gehh' 


By  introducing  notations  for  periodic  speed  and  elevation  of  the  surface  and  internal 
mode  of  oscillations,  and  by  using  equations  (3.4)  and  (3.5),  relationship  between 
current  speed  and  elevation  of  the  surface  and  internal  mode  is  found: 


u  =  u  exp(/rat) 


C  =  £  exp(; icot) ; 


and 


surface  mode: 


i— - -exp  (icot) 

co  dx 


internal  mode: 


u,  =  i  - 


g£h  dC 

co(h  +  h')  dx 


exp(icot) . 


(3-6) 


3.  2.  Progressive  surface  wave  propagating  from  deep  into  shallow  water 


Results  of  two-layer  analytical  model,  al  lowing  for  an  abrupt  depth  change  at  the  shelf 
break  and  assuming  progressive  wave  propagating  from  deep  into  shallow  water,  will 
be  considered.  Schematic  presentation  of  the  model  is  given  in  Figure  14.  Since 
model  area  has  two  parts  with  different  bottom  layer  depths  and  same  thermocline 
depth,  notation  in  equations  for  the  surface  and  internal  mode  is  slightly  changed. 
Apostrophe  which  was  denoting  influence  of  internal  mode  on  the  surface  and  bottom 
layer  will  be  replaced  by  for  negative  jc  and  for  positive  x.  The  same  notation 
will  be  used  in  the  surface  mode  equations. 
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Figure  14.  Schematic  presentation  of  idealized  topography  for  two-layer  analytical 
model  with  shelf  region  (-)  and  deep  sea  region  (+)  indicated. 


If  an  incoming  wave  with  amplitude  Ip  propagates  from  region  (+)  into  region  (-), 
abrupt  topography  change  generates  reflected  barotropic  wave  in  (+),  transmitted 
barotropic  wave  in  (-),  reflected  baroclinic  wave  in  (+)  and  transmitted  baroclinic 
wave  in  (-): 


C  =  Tp  exp[/©  +  £;x)J 

surface  mode: 

CP  =  exp[/(ft*  +  k+px) J+  R+p  exp|/(raf  -  Vpx)\ 

Cu  =  T~  exp[/(&*  +  /;ux)J 

internal  mode: 

C  =  K  exp[/(&*  -  k*x)] , 


(3.7) 


where  T~  represents  transmitted  barotropic  wave  amplitude  in  region  (-),  R+p  is 
amplitude  of  reflected  barotropic  wave  in  region  (+),  T~  denotes  transmitted 
baroclinic  wave  amplitude  in  region  (— )  and  R*u  reflected  baroclinic  wave  amplitude 
in  region  (+). 

Relationship  between  speed  and  elevation  of  the  surface  and  internal  mode  can  be 
found  on  the  basis  of  equations  (3.6): 
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surface  mode: 


u 

internal  mode: 


with: 

c~P  =  4s(h+h')  '>  cl  =  4'g(h+h*) ; 

C~  =  4g~Jt  ;  G~  =  ge  ;  cl=4&¥  ;  G+=gE-\ 

h+h  h  +  h 


(3.8) 


Wave  forms  for  the  surface  and  internal  mode  in  region  (-)  and  region  (+)  include 
four  unknown  amplitudes  which  depend  on  the  amplitude  of  incoming  surface  wave. 
Relationship  between  these  amplitudes  can  be  obtained  by  using  boundary  conditions 
at  depth  discontinuity.  At  x  =  0,  elevation  continuity  in  the  surface  and  bottom  layer, 
along  with  continuous  flow  in  the  entire  water  column  and  in  the  surface  layer,  needs 
to  be  fulfilled.  Expressions  (3.4)  and  (3.5)  are  also  taken  into  account,  showing  the 
influence  of  the  surface  mode  elevation  on  the  bottom  layer  elevation,  and  of  the 
internal  mode  speed  on  speed  in  the  surface  layer.  Boundary  conditions  are  then: 
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(3.9d) 


If  expressions  (3.7)  and  (3.8)  are  included  in  boundary  conditions  (3.9)  for  x  =  0,  four 
equations  with  five  variables  are  obtained.  In  this  way,  unknown  amplitudes  of 
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reflected  and  transmitted  waves  can  be  expressed  in  terms  of  amplitude  of  the 
incoming  surface  wave.  System  of  equations  is  as  follows: 
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(3.10b) 

(3.10c) 

(3 . 1  Od) 


From  (3.10a)  and  (3.10c),  which  are  the  same  equations  as  in  barotropic  case 
( Cerovecki ,  1992),  one  could  get: 
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(3.11) 


where  H  -h  +  h  and  H*  —  h  +  h* ,  i.e.  depths  of  shallow  and  deep  region 
respectively. 

Including  these  terms  in  equations  (3.10b)  and  (3 . 1  Od),  system  of  equations  for  R* 
and  T~  is  found.  From  (3. 1 0b): 
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Inclusion  of  term  (3.12)  in  equation  (3 . 1  Od)  gives  equation  for  R* 
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By  using  terms  for  c  ,  c+ ,  cu ,  cj,  G  and  G+ ,  the  solution  for  R,t  is  obtained,  i.e. 

the  amplitude  of  transmitted  baroclinic  wave  in  dependence  on  depths  in  the  shallow 
and  deep  region,  thermocline  depth  and  stratification  of  water  column  (e  enters  into 
equation  through  wave  speed  in  the  bottom  layer  and  via  terms  G~  and  G+ ): 
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Simplifying  this  and  by  changing  hi  -h  =  H+  -  H  one  gets: 
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so  that  the  final  expression  for  R*  is: 
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Finally,  by  using  (3.13)  in  (3.12),  the  solution  for  Tu  is: 

rh 


77=2/ 


h 


h~Jft  V/r/r  +4ini- 


1  *ihh~ 


(3.14) 
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In  this  way,  the  final  solutions  for  equation  (3.10)  are  found: 


R+p  =lp  y[JF+4rr 


JhiJJF -Jh~U  1  A fhh 


(3.15) 


The  main  factor  contributing  to  the  amplification  of  tidal  signal  on  thermocline 
(pycnocline)  because  of  the  depth  change  is: 


p  V/r//+  +4h*Jr  4e' 


which  shows  that  thermocline  oscillations  are  more  pronounced  if  thermocline  depth 
is  larger,  if  density  difference  between  the  surface  and  bottom  layer  is  smaller,  and  if 
amplitude  of  incoming  surface  tidal  wave  is  larger.  When  shallow  region  depth  is 
increased,  intensification  decreases,  while  an  increase  of  deep  region  depth  gives 
more  significant  amplification  (Figure  15).  These  results  are  in  accordance  with 
previous  analysis  for  two-layer  analytical  models  (e.g.  Rattray,  1960). 

Surface  and  interface  elevations  can  be  found  from  expressions  connecting  solutions 
for  the  surface  and  internal  mode  (equations  are  valid  both  in  the  shallow  and  deep 
regions  of  the  model,  so  they  are  quoted  in  original  notation): 


so  that  in  the  shallow  region  (-)  and  deep  region  (+),  which  are  now  being  denoted  by 
(1)  and  (2),  sea  surface  elevation  and  thermocline  elevation  are  given  by: 
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(c\ =C= tp  expK^ + *;*)J 


(Di  =  T^—  c  +  C  =  7^—  T;  exp[i(fi#  +  *;x)]+  T~  exp[/(ft#  +  k'x)] 
h+h  h+ h 

[C)2  =CP  =1P  exp[/(fi*  + *£*)]+**  exp|/(ft#-*;x)] 

(O 2  =  ,  ^  ,  +  C  +  C  =  exp[/(^  +  *^r)]+  /?;  exp[/(<tf  -  &>)]} 

/?  +  n  n  +  h 


+  Ru  exp[/^uf  -  ku  x 


The  distribution  of  the  ratio  between  thermocline  elevation  (£")  and  sea  surface 
elevation  (0  is  given  in  Figure  15,  for  x  =  0  and  range  of  depths  of  the  shallow  and 
deep  region.  In  this  example  thermocline  depth  was  set  to  /*  =  25m,  while 
stratification  coefficient  was  taken  to  be  e  —  0.00045 .  These  parameters  correspond  to 
the  characteristics  of  stratification  in  June  2001,  when  strong  diurnal  tidal  oscillations 
were  detected  at  Lastovo,  in  the  layer  between  15  and  36  m  ( Mihanovic ,  2005). 
Thermocline  oscillations  were  significantly  correlated  with  the  surface  tidal  signal  at 
Dubrovnik.  The  deeper  the  shallow  region,  the  smaller  the  ratio  between  thermocline 
elevation  and  surface  elevation,  the  deeper  the  deep  region,  the  larger  the  ratio. 
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Figure  15.  The  distribution  of  the  ratio  between  thermocline  elevation  and  sea  surface 
elevation  for  x  =  0  and  assumed  thermocline  depth  h-  25  m,  with  stratification 
coefficient  e  =  0.00045 . 
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The  distribution  of  the  ratio  between  thermocline  elevation  and  sea  surface  elevation, 
at  x  =  0,  for  different  values  of  thermocline  depth  and  stratification  coefficient,  can 
be  seen  in  Figure  16.  South  Adriatic  shelf  profile  between  Lastovo  island  and  the 
middle  part  of  the  South  Adriatic  Pit  was  analyzed.  The  least  squares  method  gives 
the  depth  of  shallow  region  H~  =  150  m,  and  the  depth  of  deep  region  H+  =  1 170  m. 
The  influence  of  thermocline  depth  increase  and  stratification  reduction  on 
thermocline  oscillations  amplification  is  evident.  When  parameters  measured  in  June 
2001  are  taken  into  account  (red  lines  in  Figure  16),  analytical  model  with  progressive 
wave  propagating  from  deep  into  shallow  region  gives  the  ratio  of  thermocline 
oscillations  and  sea  surface  oscillations  of  about  10.  The  data  from  June  2001 
demonstrate  that  the  ratios  of  maximum  filtered  values  range  between  9  and  39  (the 
envelope  of  surface  elevation  has  smaller  amplitude  than  thermocline  elevation 
envelope).  If  only  monthly  maxima  are  taken  into  consideration,  the  ratio  is  about  29. 
Harmonic  analysis  of  Kl  tidal  components  showed  that  the  amplitude  ratio  was 
approximately  21  ( Mihanovic ,  2005). 
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Figure  16.  The  distribution  of  the  ratio  between  thermocline  elevation  and  sea  surface 
elevation  for  x  =  0  and  assumed  shallow  region  depth  H~  =  1 50  m  and  deep  region 
depth  H+  =1170  m  (profile  Lastovo  -  South  Adriatic  Pit).  Red  lines  denote 
thermocline  depth  (h  =  25m)  and  stratification  (£  =  0.00045)  as  measured  in  June 
2001  {Mihanovic,  2005). 
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These  results  imply  that  the  model  defined  in  this  way  can  reproduce  relationship 
between  thennocline  elevation  and  surface  elevation,  despite  the  model  limitations 
(two-layer  analytical  model,  abrupt  topography  change  at  the  shelf  break).  However, 
in  progressive  wave  consideration,  surface  oscillations  and  current  oscillations  are  in 
phase,  which  does  not  comply  with  K1  tidal  component  characteristics  in  the  Adriatic. 
Measurements  from  July  1990  and  August  1991  showed  that  this  relationship  is  closer 
to  standing  surface  wave  character,  since  current  oscillations  precede  surface 
oscillations  by  about  90°  ( Mihanovic ,  2005).  Therefore,  the  model  results  with  the 
standing  surface  wave  will  be  considered  in  the  next  section. 


3.3.  Standing  surface  wave 

Harmonic  analysis  of  sea  currents  in  a  wider  Lastovo  area  showed  that  maxima  of 
barotropic  flow  precede  maxima  of  sea  level  by  about  90°.  Analytical  model  with 
progressive  surface  wave  cannot  reproduce  this  phase  relationship,  since  for 
progressive  wave  the  flow  and  sea  level  oscillations  are  in  phase.  Consequently,  two- 
layer  analytical  model  with  assumption  of  a  standing  surface  wave  was  also 
developed.  As  in  previous  section,  elevation  continuity  condition  in  the  surface  and 
bottom  layer,  as  well  as  flow  continuity  in  the  entire  water  column  and  separately  in 
the  surface  and  bottom  layer,  needs  to  be  fulfilled. 

Generally,  system  of  equations  for  the  surface  (3.4)  and  internal  mode  (3.5)  can  be 
written  in  transformed  form,  which  uses  surface  mode  flow  and  internal  mode  flow  in 
the  surface  layer,  instead  of  speed  of  the  surface  and  internal  mode.  The  notation  is 
original,  i.e.  h  is  depth  of  the  surface  layer,  h'  is  depth  of  the  bottom  layer,  u  and  £ 
represent  speed  and  elevation  of  the  surface  layer,  while  u  and  represent  speed 
and  elevation  of  the  bottom  layer.  If  speed  and  elevation  of  the  surface  mode  are 
denoted  by  u  and  C,  ,  and  speed  and  elevation  of  the  internal  mode  by  uu  and  £u , 

new  variables  can  be  defined,  which  will  be  used  in  transformed  system  of  equations. 
These  new  variables  are  as  follows: 

■  U p  =(h  +  h')u  -  surface  mode  flow, 

■  C p  =  £ P  ~  surface  mode  elevation  (equals  surface  elevation), 

■  U u  =  huu  -  internal  mode  flow  in  the  surface  layer,  and 

■  Cu  =  ~CU  '  transformed  internal  mode  elevation  (enables  similar  equations 
for  both  modes). 


Equations  for  the  surface  and  internal  mode  now  become: 


surface  mode: 


+  c ; 


K 

dx 


=  0 


3 1  dx 


=  0 


where  c2  =g(h  +  h'). 


(3.16) 
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internal  mode: 


,  2  hh' 

where  cu  =  g£ - . 

h  +  h' 


(3.17) 


dU„ 


+  cl^  =  0 


dt  “  dx 

*£u  ,  at/. 


-  +  - 

dt  dx 


=  0 


Boundary  conditions  at  depth  discontinuity  are  equal  elevations  and  flows  at  x  =  0, 
giving: 

surface  mode:  Up  and  £p  are  preserved  at  depth  discontinuity,  (3.18a) 


internal  mode: 


—  £u  H - is  preserved  at  depth  discontinuity, 

h  +  h 

(3.18b) 

h 

- U  +  U  is  preserved  at  depth  discontinuity. 

h  +  h'  p 


The  relationship  with  the  original  variables  (speed  and  elevation  of  the  surface  layer, 
and  speed  and  elevation  of  the  bottom  layer)  is  defined  by: 


u=  1  ,UD+-  Uu 
h  +  h  h 

(3.19a) 

VTNf 

II 

(3.19b) 

u'  =  1  U  -  ]  U 

h  +  h'  p  h'  u 

(3.19c) 

h' 

r=  n  r  ~r . 

^  h+h,bp  " 

(3. 1 9d) 

In  the  following  equations,  all  variables  to  the  left  of  the  origin  (depth  discontinuity) 
will  be  denoted  by  and  those  to  the  right  by  so  that  the  apostrophe  mark  is 
temporarily  lost.  Therefore,  depths  of  the  bottom  layer  on  the  left  and  on  the  right  side 
of  topography  jump  become,  as  in  section  3.2,  h~  and  h+ .  Additional  abbreviation  is 
introduced: 


b  h(h'-h~) 

\h  +  tfjh  +  *-)' 
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Generally,  the  solution  will  be: 


X(x,t) 


~Vp{x,tJ 

Cp{xj) 
uu  (x,/) 

=  X{x]e,aM  = 

f,W 

CM 

.CMO. 

_CM_ 

meaning  that  the  solution  has  a  defined  angular  frequency  co.  The  solutions  of 
equations  (3.16)  and  (3.17)  can  be  written  separately  for  the  left  and  right  side.  For  the 
left  side: 


u;(xj 

CO 

cos — X 

i  .  co 

.  _  .  CO 
—  ic  sin  —  x 

c; 

CO 

0 

0 

[mmi 

CM 

- sin — x 

c  c 

p  p 

cos — X 

0 

0 

CM 

u;(x) 

b  co 

— cos — X 

ib  _  .  co 
— c,  sin — x 

CO 

cos  —  X 

.  CO 

—  ic  sm  —  x 

u[:r)(  o) 

CM. 

2 

ib  .  co 

- sin — x 

2c  c 

L  M  U 

2 

b  co 

—  cos  —  X 

2  C 

c; 

i  .  co 

- sin  —  x 

c,  c 

u  u 

U  — 

Cu 

CO 

cos — X 

C~u 

CM(0). 

where  U  ( 0),  CM  represent  values  of  flow  and  elevation  of  the  surface  mode  at 

x  =  0  (preserved  on  depth  discontinuity),  while  0),  CHo)  represent  average 

values  of  corresponding  parameters  for  the  internal  mode  on  both  sides  of 
discontinuity  (not  preserved).  These  four  complex  numbers  completely  determine  the 
solution  of  the  form  X {x^"* . 

For  the  right  side  the  solution  is: 


X+{x)  = 


u;(x) 

CM 

iCto 

.CM 


a ) 

cos — x 

r  + 
c  P 

i  .  co 

- sin  — x 

c+  c+ 

p  p 

b  co 

—cos — X 
2  cl 

ib  .  co 
- sin  —  x 

2c:  c: 


.  +  .  co 
—  icn  sin  —  x 

'  c; 

co 

cos  —  X 

c 
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2 
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—cos  —  X 
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co 

cos  —  X 

c 

/  .  CO 

— sin  —  x 

c+  c+ 


.  +  .  co 
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u 

co 
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CM 
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This  solution  is  proportional  to  e'°* ,  i.e.  it  is  related  to  frequency  co.  Frequency  -co 
has  complex  conjugate  solution,  i.e.  X~  and  X+ .  Therefore,  respective  constants 
U  (o),  ...,  C'T,(0)  f°rrn  complex  conjugate  pairs  and  by  combining  both  solutions 


(3.20) 


(3.21) 
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one  can  find  real  solution  characterized  by  frequency  ca.  Thermocline  elevation 
(3. 1 9d)  contains  reduced  surface  mode  elevation  and  internal  mode  elevation,  which 
include  four  terms  in  the  last  row  of  (3.20)  and  (3.21).  The  amplification  of 
thermocline  elevation  arises  form  terms  which  multiply  U  (6)  and  ^'^(0) .  It  is 

obvious  that  term  connected  with  U  { 0)  is  divided  by  the  speed  of  internal  waves  c*, 
which  is  significantly  smaller  than  the  speed  of  surface  waves,  contributing  to 
amplification  on  thermocline.  Intensification  related  to  0)  is  not  seen  directly, 

because  it  comes  from  the  constant  already  having  amplified  thermocline  elevation.  It 
becomes  more  comprehensible  when  boundary  conditions  at  the  ends  of  the  basin  are 
included,  thus  enabling  the  boundary  problem  to  be  solved.  On  the  other  hand,  this 
term  is  essential  for  reproducing  propagating  character  of  internal  waves.  Therefore, 
two  different  cases  of  solution  with  boundary  conditions  taken  into  account  will  be 
presented  in  this  report.  First,  the  problem  will  be  solved  with  classical  boundary 
conditions  on  the  left  and  right  end  of  the  basin  for  the  surface  and  internal  mode  (at 
the  closed  boundary  no  normal  flow  condition,  and  at  the  open  boundary  zero 
elevation  condition).  After  that,  the  analysis  of  a  problem  in  which  surface  mode  is 
subjected  to  classical  boundary  conditions  and  internal  mode  to  radiation  boundary 
condition  will  be  offered  (so  that  the  boundaries  of  the  basin  do  not  affect  wave 
disturbance,  allowing  complete  energy  transmission).  Schematic  presentation  of  a 
basin  that  will  represent  the  Adriatic  is  given  in  Figure  17.  The  left  end  of  the  basin  at 
x  =  L (  <  0  is  closed,  whilst  the  right  end  at  x  =  L2  >  0  is  open,  i.e.  it  borders  with  an 
ocean  of  infinite  depth. 
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Figure  17.  Schematic  presentation  of  an  idealized  topography  for  two-layer  analytical 
model  with  shelf  region  (-)  and  deep  sea  region  (+).  In  the  shelf  region  (-)  the  basin 
is  closed,  and  in  the  deep  sea  region  (+)  it  borders  with  an  ocean  of  infinite  depth. 
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In  order  to  implement  the  analytical  model  to  the  Adriatic,  it  is  necessary  to  determine 
parameters  that  approximate  the  Adriatic  Sea  with  a  two-dimensional  basin  closed  at 
left  end,  and  opened  at  the  right  end,  with  rapid  depth  change  at  the  shelf  break. 
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Parameters  which  are  included  in  the  model  are  depths  of  the  shallow  and  deep  region 
(H~  and  H+ ),  together  with  distance  of  the  closed  (Z,)  and  open  end  (Z2 )  from  the 
origin.  These  parameters  were  resolved  in  three  ways.  The  first  approximation  that 
was  used  is  the  Adriatic  transversal  cross  section  approximation,  with  sections  that  are 
20.5  km  apart  ( Sterneck ,  1919),  which  enabled  calculation  of  average  depths  of  these 
sections  along  the  Adriatic.  Next,  these  values  separated  by  20.5  km  were  interpolated 
every  0.5  km  enabling  more  precise  discretization  and  better  resolving  of  depth 
discontinuity  position.  The  third  Adriatic  approximation  is  based  on  longitudinal  cross 
section  along  central  Adriatic  line,  also  interpolated  every  0.5  km. 

Three  different  parameterizations  of  the  Adriatic  are  shown  in  Figure  18: 

-  parameterization  according  to  Sterneck,  spacing  equals  20.5  km, 

-  parameterization  according  to  Sterneck,  spacing  equals  0.5  km, 

-  central  longitudinal  Adriatic  cross  section,  spacing  equals  0.5  km. 
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Figure  18.  The  Adriatic  Sea  approximated  by  a  two-dimensional  basin  with  three 
different  parameterizations.  Green  line  denotes  geographical  position  of  Lastovo  in 
idealized  basin  representing  the  Adriatic. 


The  exact  position  of  the  depth  discontinuity,  together  with  depths  of  the  shallow  and 
deep  region,  was  resolved  by  applying  the  least  squares  method  on  depth  values 
distributed  along  the  Adriatic.  Lastovo  island  is  located  around  490  km  from  the 
closed  end  of  the  basin  (represented  by  green  line  in  Figure  18).  Once  the  coordinate 
of  the  discontinuity  is  determined,  origin  is  positioned  at  this  coordinate,  and  model 
parameters  are  defined.  They  are  given  in  Table  2,  for  each  parameterization  in 
respective  coordinate  system.  Depths  are  smaller  for  Sterneck’s  approximations 
because  they  comprise  average  depth  of  each  transversal  transect,  whereas 
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longitudinal  Adriatic  cross  section  includes  maximum  depths  on  transversal  cross 
sections. 

In  the  three  cases  we  have  used  thermocline  depth  and  stratification  coefficient  as 
calculated  for  June  2001  and  Lastovo  (h  =  25  m,  e  =  0.00045  ,  Mihanovic,  2005). 


Table  2.  Parameters  required  for  the  two-dimensional  analytical  model  including  a 
standing  surface  wave  and  three  different  approximations  of  the  Adriatic. 


Sterneck 

20.5  km 

Sterneck 

0.5  km 

Central 

longitudinal  cross 
section  of  the 
Adriatic 

A 

-533  km 

-540.5  km 

-514  km 

^2 

267  km 

259.5  km 

286  km 

H~ 

76.5  m 

77.6  m 

111.9m 

H+ 

491.7  m 

489.9  m 

1014.5  m 

Lastovo 

-43  kin 

-50.5  km 

-24  km 

3.  3.  1.  Classical  boundary  conditions  at  the  closed  and  open  end  of  the  basin  for  the 
surface  and  internal  mode 

Classical  boundary  conditions  for  the  surface  and  internal  mode  are  basic  assumptions 
of  the  model,  from  which  relationship  between  constants  U  ( 0),  ...,  Ct'K 0)  can  be 

found,  i.e.  connection  between  surface  oscillations  and  thermocline  response  may  be 
determined.  At  the  closed  end  of  the  basin,  corresponding  flows  are  zero,  as  are 
elevations  at  the  open  end: 

closed  (left)  end  of  the  basin:  Up(Ll)  =  0,  Uu(Li)  =  0,  (3.22a) 

open  (right)  end  of  the  basin:  £u(L2)  =  0.  (3.22b) 


It  is  necessary  to  resolve  the  frequency  a),  together  with  constants  (/  (0),  f p (0) , 
U^( 0)  and  C^(0),  so  that  solutions  (3.20)  and  (3.21)  comply  with  boundary 
conditions  (3.22).  In  the  following  equations  these  constants  will  be  denoted  by  ylf 
y2 ,  y3  and  y4 .  Non-trivial  solutions  of  system  (3.23)  need  to  be  found: 
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Non-trivial  solutions  exist  for  those  co  for  which  determinant  of  the  system  is  zero. 


i.e.: 
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If  the  term  in  the  first  bracket  equals  zero,  natural  frequencies  of  surface  modes  are 
found,  whereas  the  second  bracket  equaling  zero  gives  internal  modes  frequencies. 
Natural  frequencies  of  surface  modes  essentially  represent  characteristic  frequencies 
of  fundamental  mode  and  higher  seiche  modes  for  basin  defined  by  depths  and 
distances  from  the  closed  and  open  basin  boundary.  For  these  frequencies,  one  needs 
to  solve  system  (3.23).  For  internal  modes  y,  =  y2  = 0 ,  so  that  they  are  manifested 
only  through  and  Uu .  For  surface  modes,  all  four  parameters  are  non-zero, 
meaning  that  surface  modes  influence  all  four  variables,  i.e.  surface  and  internal 
elevation  and  speed.  It  is  important  to  emphasize  that  y}  and  y4  depend  on  yl  and 

y2  in  a  way  which  depends  on  given  boundary  conditions. 

If  the  lower  right  2x2  matrix  in  (3.23)  is  denoted  by  Du ,  i.e. 
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the  relationship  between  y3  and  y4  as  well  as  y]  and  y2  could  be  expressed  as: 
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or: 
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Also,  from  equation  (3.23),  for  <a  which  cancels  out  the  left  bracket  in  (3.24),  one  can 
get: 
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which  together  with  previous  equation  gives: 
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By  substituting  solutions  yt,  y2>  y2  and  y$  with  £/,(o),  ^(o),  z/,(;r)(o)  and 
C>(0)  in  (3.20),  the  solution  for  the  left  side  of  the  basin  in  matrix  form  is  found 
(the  solution  for  the  right  side  is  analogous,  according  to  (3.21)): 


(3.25) 
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By  arranging  elevation  terms  on  the  left  side  and  on  the  right  side  of  the  origin,  final 
solutions  that  enable  comparison  of  oscillations  in  the  surface  and  bottom  layer  are 
obtained: 
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Surface  layer  elevation  and  bottom  layer  elevation  are  obtained  from  solutions  for  the 
surface  and  internal  mode  according  to  (3.19). 

General  solutions  (3.20)  and  (3.21)  demand  designation  of  flow  and  elevation  for  the 
surface  mode  (f/p(0),  <^(o)),  as  well  as  average  values  of  flow  and  elevation  for  the 

internal  mode  ( U^'r> (0) ,  0))  at  depth  discontinuity.  If  the  average  values  for  the 

internal  mode  are  disregarded,  i.e.  if  it  is  supposed: 
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amplification  for  is  obtained  only  through  the  sine  term  (factor  cpjcu  multiplies 

sine  in  the  lower  left  element  of  (3.20)  and  (3.21)).  This  ignores  the  influence  of 
surface  mode  on  internal  mode,  which  emerges  clearly  when  solving  the  boundary 
problem,  where  surface  mode  values  (_y,  and  y2 )  determine  values  for  internal  mode 
(y2  and  y4 ).  This  is  explicitly  clear  in  equation  (3.25).  Factor  c pjcu  is  evident  in  the 

term  for  y4 ,  and  comes  from  yx  (representing  the  speed  at  depth  discontinuity). 
Hence,  the  speed  at  topography  discontinuity  in  the  internal  mode  is  responsible  for 
amplification  of  thermocline  oscillations.  Factor  cpjcu  in  the  final  solution 

represents  this  intensification,  and  it  can  be  seen  in  sine  term  (through  yx )  and  in 
cosine  term  (through  y4 ),  which  is  in  accordance  with  progressive  character  of 
internal  waves. 

The  results  of  analytical  model  for  the  surface  and  bottom  layer,  with  classical 
boundary  conditions  at  the  closed  and  open  end  of  the  basin,  are  given  in  Figures  19  - 
21.  Three  parameterizations  of  the  Adriatic  were  used,  corresponding  to  a  two- 
dimensional  basin  with  a  depth  discontinuity  (Figure  18  and  Table  2).  Geographical 
position  of  Lastovo  island  is  marked  by  a  green  line.  Since  the  solutions  given  in 
(3.27)  are  defined  up  to  an  arbitrary  constant,  results  with  the  supposed  surface 
amplitude  of  1  m  at  the  closed  end  of  the  basin  are  presented. 
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Figure  19.  The  results  of  the  two-layer  analytical  model  for  the  surface  and  bottom 
layer  with  classical  boundary  conditions  at  the  closed  and  open  end  of  the  basin,  and 
parameterization  of  the  Adriatic  according  to  Sterneck,  with  spacing  of  20.5  km.  The 
position  of  Lastovo  is  marked  by  a  green  line. 
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Figure  20.  The  results  of  the  two-layer  analytical  model  for  the  surface  and  bottom 
layer  with  classical  boundary  conditions  at  the  closed  and  open  end  of  the  basin,  and 
parameterization  of  the  Adriatic  according  to  Sterneck,  with  spacing  of  0.5  km.  The 
position  of  Lastovo  is  marked  by  a  green  line. 


-500  -400  -300  -200  -100  0  100  200 


Figure  21.  The  results  of  the  two-layer  analytical  model  for  the  surface  and  bottom 
layer  with  classical  boundary  conditions  at  the  closed  and  open  end  of  the  basin,  and 
parameterization  of  the  Adriatic  along  central  longitudinal  cross  section,  with  spacing 
of  0.5  km.  The  position  of  Lastovo  is  marked  by  a  green  line. 
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By  solving  equations  (3.24)  natural  frequencies  of  the  surface  modes  are  calculated. 
They  represent  characteristic  frequencies  of  fundamental  mode  and  higher  seiche 
modes  in  the  basin  defined  by  depths  of  shallow  and  deep  region,  as  well  as  by 
distances  of  closed  and  open  boundary  from  the  origin.  Parameterization  of  the 
Adriatic  according  to  Sterneck  gives  period  of  the  fundamental  mode  of  about  23.34  h 
(for  spacing  of  20.5  km)  and  23.46  h  (for  spacing  of  0.5  km).  When  central 
longitudinal  Adriatic  approximation  is  used,  the  fundamental  mode  period  is  about 
18.32  h.  These  estimates  are  close  to  measured  values  of  fundamental  seiche  period  in 
the  Adriatic,  which  lie  between  21.4  and  22.4  h  (Buljan  and  Zore-Armanda ,  1976; 
Mosetti  and  Purga ,  1983),  taking  into  account  that  these  are  approximations  of  the 
Adriatic  by  a  simple  two-dimensional  analytical  model.  Estimates  calculated  from 
Sterneck’s  parameterizations  are  closer  to  cited  periods,  since  depth  values  represent 
average  values  at  a  transversal  cross  section.  In  the  central  longitudinal  Adriatic  cross 
section,  only  maximum  depths  are  used,  giving  smaller  fundamental  mode  period  for 
this  approximation. 

All  three  parameterizations  give  amplification  of  thermocline  oscillations  with  respect 
to  surface  oscillations.  However,  the  ratio  of  amplitudes  as  well  as  phase  relationship 
sensitively  depend  on  selection  of  depth  parameters  (shal  low  and  deep  region)  and  of 
distances  of  the  closed  and  open  boundary  from  the  origin.  Since  solutions  for  the 
surface  and  bottom  layer  are  in  the  form  of  standing  waves,  and  since  lengths  of 
internal  waves  are  much  smaller  than  those  of  surface  waves,  oscillations  can  be 
either  in  phase  or  opposed  in  phase.  In  the  case  of  Stemeck’s  parameterization  with 
spacing  of  20.5  km  (Figure  19),  significant  amplification  is  detected  at  Lastovo 
position  (about  100  times).  Yet,  oscillations  are  opposed  in  phase,  which  is  not  in 
accordance  with  observed  “in  phase”  relationship  ( Mihanovic ,  2005).  A  slight  change 
of  parameters  for  Sterneck’s  0.5  km  parameterization  gives  completely  different 
results  (Figure  20).  The  amplification  of  oscillations  is  more  pronounced  in  deep 
region  of  the  basin,  and  Lastovo  is  positioned  almost  at  the  nodal  point  of  internal 
standing  waves.  This  results  in  insignificant  thermocline  oscillations  at  this  location 
(the  ratio  between  thermocline  elevation  amplitude  and  surface  elevation  amplitude  is 
about  0.2).  Finally,  parameterization  of  the  Adriatic  based  on  central  longitudinal 
cross  section  also  gives  significant  amplification  in  both  shallow  and  deep  region,  but 
again  it  shows  sensitive  dependence  of  the  model  on  values  of  calculated  parameters 
(Figure  21).  More  precisely,  at  Lastovo  the  thermocline  and  surface  oscillations  are  in 
phase,  and  amplification  ratio  is  about  190. 

Considering  that  in  the  surface  layer  standing  wave  is  presumed,  and  boundary 
conditions  give  maximum  surface  elevation  at  the  closed  boundary,  surface 
oscillations  have  relatively  small  amplitude  close  to  the  origin  (depth  discontinuity) 
when  compared  with  the  oscillations  near  the  closed  boundary.  For  this  reason, 
different  basin  parameterizations  give  unrealistic  amplifications  of  thermocline 
oscillations.  The  ratios  of  thermocline  and  surface  oscillations  as  well  as  phase 
differences  between  them  depend  on  the  bottom  layer  results  and  on  the  Lastovo 
position  with  respect  to  nodal  points  of  internal  standing  waves.  Therefore,  it  can  be 
concluded  that  this  analytical  model  explains  to  some  extent  physical  origin  of 
amplification  of  thermocline  oscillations  at  depth  discontinuity  (velocity  change,  or 
better  to  say,  barotropic  flow  change).  However,  it  fails  to  reproduce  phase  and 
amplitude  relationships  as  recorded  in  June  2001. 


41 


3.  3.  2.  Classical  boundary  conditions  for  the  surface  mode  and  radiation  boundary 
condition  for  the  internal  mode  at  the  closed  and  open  end  of  the  basin 

In  previous  section  classical  boundary  conditions  for  the  surface  and  internal  mode 
were  presumed.  This  section  will  deal  with  radiation  boundary  conditions  for  internal 
mode: 


closed  (left)  end  of  the  basin:  Up  \LX )  =  0 , 

open  (right)  end  of  the  basin:  C  P (A ) =  0 » 
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General  solution  of  system  of  equations  (3.16)  and  (3.17)  can  be  written  as: 
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where  following  matrices  were  introduced: 
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Unknown  coefficients  ylt  y2 ,  y3  and  y4  can  be  obtained  by  applying  boundary 
conditions,  which  leads  to  the  system  of  equations  analogous  to  system  (3.23): 
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were  following  matrices  were  introduced: 
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Natural  frequencies  are  solutions  of  the  equation: 

det(Dp)det(D.)=0, 


(3.32) 


where  left  factor  corresponds  to  surface  modes,  the  right  one  to  internal  modes.  After 
obtaining  frequency  of  the  fundamental  mode,  from  (3.23)  it  follows  as  before: 
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along  with: 
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It  is  not  necessary  to  explicitly  calculate  coefficients  Yu,  because  previous  equation 
can  be  included  into  solutions  (3.29)  and  (3.30): 
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Finally,  the  terms  which  need  to  be  written  in  complete  form,  introduced  in  (3.33), 
are: 


- Y  +Y  =  ~bD 

2  P  u  u 


'0  0" 

D  Y„, 

L 

-Y „+Y  =b  ZT1 

'l 

o' 

0  1 

u  p  “ 

2  II  M 

0 

0 

d.  r. 


(3.34) 


The  part  of  the  solution  corresponding  to  surface  mode  is  the  same  as  in  previous 
section,  while  the  solution  for  internal  mode  states: 
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where  determinant  det (Du )  is  given  by: 
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The  general  solution  for  free  oscillations  with  angular  frequency  co,  as  a  function  of 
time  and  space  (real  form)  and  up  to  arbitrary  amplitude^  and  phase  cp,  is: 
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This  model  was  also  combined  with  three  parameterizations  of  the  Adriatic,  as 
described  in  the  previous  section.  The  most  important  difference  arising  from  the 
model  with  radiation  boundary  conditions  for  internal  mode  is  the  progressive 
character  of  waves  on  thermocline.  Amplification  of  tidal  signal  occurs  due  to  the 
same  reasons  as  in  the  previous  case,  i.e.  due  to  the  ratio  between  speed  of  the  surface 
and  speed  of  the  internal  waves.  However,  this  model  is  also  sensitively  dependent  on 
model  parameters  -  depths  of  shallow  and  deep  region  and  discontinuity  coordinate. 
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The  coordinate  corresponding  to  the  position  of  Lastovo  is  also  very  important  for  the 
results  of  the  model,  having  in  mind  that  surface  wave  is  standing  wave,  and  internal 
wave  has  a  progressive  character,  with  phase  depending  upon  position.  For  this 
reason,  matching  phases  between  surface  oscillations  and  thermocline  oscillations  will 
occur  only  at  a  certain  point,  its  position  depending  on  the  speed  of  internal  waves. 
The  best  results  were  found  for  Sterneck’s  approximation,  with  spacing  equal  to  20.5 
km,  and  these  results  are  presented  in  Figures  22  -  23. 

Figure  22  shows  space-time  distribution  of  surface  elevation  for  a  distance  extending 
100  km  from  the  origin  in  both  directions,  according  to  equations  (3.27a)  and  (3.27c). 
The  time  interval  of  48  h  was  analyzed,  with  natural  frequency  of  the  fundamental 
mode  related  to  above  mentioned  Adriatic  parameterization.  Standing  wave  can  be 
noticed,  with  larger  amplitudes  obtained  to  the  left  of  the  origin,  which  corresponds 
well  with  the  results  in  the  previous  section  (Figures  19  to  21).  Space-time 
distribution  of  thermocline  elevation  (Figure  23)  was  calculated  by  including  terms 
(3.27a),  (3.27c)  and  (3.35)  into  equation  (3.19d),  and  by  taking  into  account 
fundamental  frequency  for  Stemeck’s  approximation,  with  spacing  of  20.5  km. 
Thermocline  oscillations  are  amplified  with  respect  to  surface  oscillations  and  internal 
wave  is  progressive  on  both  sides  of  the  origin.  When  Lastovo  position  is  analyzed 
( x  =  -43  km),  it  can  be  seen  that  surface  and  thermocline  oscillations  are  in  phase  at 
that  position,  with  amplification  factor  of  about  50  (Figure  24).  These  results  accord 
well  with  temperature  data  measured  at  Lastovo  in  June  2001  ( Mihanovic ,  2005). 
Still,  small  change  of  depths,  coordinate  of  discontinuity  and  L.astovo  position  can 
significantly  change  the  phase  and  amplification  for  such  an  approximation. 


100  BO  60  40  JO  O  20  40  60  90  100 


Distance  (km) 


OS  OS  04  02  0  02  04  OS  OB 


Surface  elevation  (m) 

Figure  22.  Space-time  distribution  of  surface  elevation  for  a  distance  extending  100 
km  from  the  origin  in  both  directions  and  time  interval  of  48  h  (parameterization  of 
the  Adriatic  according  to  Stemeck,  with  spacing  of  20.5  km). 
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Figure  23.  Space-time  distribution  of  thermocline  elevation  for  a  distance  extending 
100  km  from  the  origin  in  both  directions  and  time  interval  of  48  h  (parameterization 
of  the  Adriatic  according  to  Sterneck,  with  spacing  of  20.5  km). 

3.  4.  Discussion 

Two-layer  analytical  model  represents  a  starting  point  for  theoretical  study  of  the 
generation  of  internal  tides.  Its  application  to  the  Adriatic  enables  one  to  acquire 
fundamental  understanding  of  physical  processes  important  for  amplification  of 
thermocline  oscillations.  It  also  paves  the  way  to  future  examination  of  internal  tides 
with  more  advanced  theoretical  models,  which  would  encompass  cross-basin  changes 
of  barotropic  flow  and  realistic  density  distribution.  Another  important  aspect  is  its 
relevance  for  development  of  numerical  models. 

This  report  examines  the  results  of  two  analytical  models  based  on  two-layer  fluid 
approximation,  with  a  rapid  depth  change  at  the  shelf  break.  The  two  models  originate 
in  the  same  initial  equations  and  utilize  the  same  boundary  conditions  at  depth 
discontinuity.  Boundary  conditions  include  elevation  continuity  in  the  surface  and 
bottom  layer,  as  well  as  flow  continuity  in  the  entire  water  column  and  in  each  layer. 
The  main  difference  between  analyzed  models  arises  from  assumption  about  the 
nature  of  the  surface  wave,  which  represents  barotropic  tide  of  diurnal  period. 
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Figure  24.  Time  series  of  surface  elevation  multiplied  by  50  (red  line)  and  time  series 
of  thermocline  elevation  (blue  line)  at  a  distance  of  34  to  50  km  from  the  origin 
towards  the  closed  end  of  the  basin  (parameterization  of  the  Adriatic  according  to 
Sterneck,  with  spacing  of  20.5  km).  Graph  representing  approximate  position  of 
Lastovo  (x-  -43  km)  is  denoted  by  an  arrow. 
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The  first  model  presumes  that  K1  tide  in  the  Adriatic  is  progressive  barotropic  wave 
propagating  from  deep  into  shallow  water.  It  is  then  partly  reflected  and  partly 
transmitted  at  topography  discontinuity,  thus  generating  reflected  and  transmitted 
internal  waves.  The  results  of  the  model  for  depth  parameters  in  the  shallow  and  deep 
region  valid  for  the  Adriatic,  and  for  thermocline  depth  and  stratification  observed  in 
June  2001,  are  close  to  the  results  determined  from  the  surface  and  thermocline 
elevations  measured  in  the  same  interval.  They  are  also  comparable  to  the  ratio  of 
amplitudes  calculated  by  harmonic  analysis  of  surface  elevation  at  Dubrovnik  and 
thermocline  elevation  at  Lastovo  for  June  2001 .  Amplitude  of  thermocline  oscillations 
is  proportional  to  surface  oscillations  and  thermocline  depth,  and  inversely 
proportional  to  density  difference  between  the  surface  and  bottom  layer.  By 
increasing  the  depth  of  shallow  region  amplification  of  thermocline  oscillations  is 
decreased,  whereas  an  increase  of  deep  region  depth  has  an  opposite  effect.  These 
results  are  in  agreement  with  previous  findings  from  similar  two-layer  analytical 
models.  However,  the  progressive  character  of  surface  wave  presumed  in  this  models 
means  that  elevation  and  velocity  oscillations  are  in  phase,  which  is  not  the  case  in  a 
wider  Lastovo  area,  or  in  accordance  with  standing  wave  character  of  diurnal  tide  in 
the  Adriatic.  Harmonic  analysis  of  sea  currents  in  the  Lastovo  channel  (July  1990) 
shows  that  diurnal  tidal  currents  lead  elevation  oscillations  by  about  90°.  This  implies 
that  surface  standing  wave  is  more  suitable  for  describing  the  characteristics  of  K1 
tide  than  is  surface  progressive  wave. 

Therefore,  another  analytical  model  is  developed,  which  also  specifies  conditions  at 
the  basin  boundaries,  thus  enabling  examination  of  surface  standing  wave  and  its 
influence  on  thermocline  oscillations.  In  the  first  case,  classical  boundary  conditions 
for  the  surface  and  internal  mode  are  presumed  at  the  closed  and  open  end  of  the 
basin.  Model  results  demonstrate  sensitive  dependence  on  the  values  of  parameters 
defining  the  Adriatic  (depths  of  shallow  and  deep  region,  distance  of  the  boundaries 
from  the  depth  discontinuity).  Since  the  solutions  for  the  surface  and  boundary  layer 
are  standing  waves,  and  lengths  of  internal  waves  lie  between  20  and  40  km,  at  some 
points  along  the  Adriatic  transect  thermocline  and  surface  oscillations  are  in  phase 
and  at  other  points  they  are  opposed  in  phase.  Depending  on  parameterization, 
Lastovo  could  be  in  either  of  these  two  areas,  and  moreover  it  could  be  positioned 
close  to  the  nodal  point.  If  it  were  at  the  nodal  point,  no  significant  amplification  of 
thermocline  variability  would  be  recorded.  When  radiation  boundary  conditions  are 
imposed  for  internal  mode,  thermocline  wave  has  progressive  character.  This  is  closer 
to  basic  ideas  about  generation  of  internal  tides,  which  state  that  barotropic  flow 
changes  on  topography  discontinuity  generate  internal  waves  that  propagate  towards 
shallow  and  deep  regions  of  the  analyzed  basin  ( Baines ,  1982).  However,  the 
sensitive  dependence  of  the  model  on  small  changes  in  parameters  defining  the 
Adriatic,  especially  on  depths  of  the  shallow  and  deep  region,  is  still  present.  The 
phase  locking  of  the  surface  and  thermocline  elevation  as  observed  at  Lastovo  is 
obtained  at  individual  points  only,  depending  on  their  distance  from  the  origin  (depth 
discontinuity)  and  the  speed  of  internal  wave.  The  best  results  are  obtained  for  the 
coarsest  Adriatic  approximation,  because  Lastovo  position  in  that  case  matches 
position  of  the  point  at  which  internal  wave  in  the  shallow  region  has  the  same  phase 
as  the  surface  standing  wave. 
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Implementation  of  the  models  described  enables  one  to  get  some  insight  into  physical 
processes  responsible  for  generation  of  internal  tides  in  the  Adriatic.  It  also  allows 
comparison  of  model  results  with  the  data.  Limitations  of  models  are  detected  and 
discussed,  being  related  to  the  character  of  surface  wave  in  the  first  model,  to  the 
Adriatic  parameterizations  and  boundary  conditions  in  the  second  model. 
Nevertheless,  it  is  hoped  that  the  present  findings  could  provide  a  starting  point  for  a 
further  investigation  of  the  Adriatic  internal  tides. 
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4.  SUMMARY  AND  CONCLUSION 


Analysis  of  historical  data  has  shown  that  in  the  middle  Adriatic  area  seasonal 
variations  of  main  atmospheric  parameters  such  as  solar  radiation,  wind,  air 
temperature  and  precipitation  cause  pronounced  seasonality  of  heat  and  water  fluxes 
between  the  atmosphere  and  the  sea.  Strong  cyclonic  activity,  especially  in  the  colder 
part  of  the  year,  causes  different  weather  types  to  occur,  with  frequent  exchange  of 
Sirocco  and  Bora  events.  During  summer,  this  area  is  under  the  influence  of 
subtropical  high-pressure  belt,  and  consequently  warm  and  dry  weather  conditions 
prevail  with  well-developed  sea-land  breeze  system.  Seasonal  cycle  of  sea 
temperature  at  various  depths  of  the  Split-Gargano  transect  points  to  isothermal 
conditions  in  late  autumn  and  winter  and  rather  small  changes  at  100  meters  where  the 
temperature  does  not  fall  below  12°C.  In  deeper  layers  annual  maximum  is  lagging 
the  surface  one,  whereas  minimum  at  all  depths  occurs  within  the  same  month.  Sea 
temperature  variations  in  a  particular  month  are  maximal  in  the  surface  layer,  with  the 
variability  decreasing  as  the  depth  increases.  Seasonal  cycle  of  salinity  shows  that  the 
salinity  increases  with  depth  throughout  the  year.  In  the  surface  layer,  down  to  10  meter 
depth,  there  are  two  minima  of  salinity,  in  May  and  July.  At  a  20  meter  depth  the  second 
minimum  disappears  whereas  the  summer-autumn  maximum  is  well  visible.  Large 
salinity  fluctuations,  maximal  in  the  surface  layer,  are  apparent  within  a  month. 

A  simple  two-layer  analytical  model  with  step-like  topography  has  been  built  as  a  first 
step  in  theoretical  analysis  of  internal  tides.  Basic  assumptions  underlying  such  a 
model  are  fulfilled  in  the  vicinity  of  the  Palagruza  Sill.  The  model  aims  at 
understanding  origin  of  internal  tides  and  it  is  based  on  separate  equations  for  the 
barotropic  and  baroclinic  mode  of  oscillations  in  the  two-layer  fluid  contained  in 
rectangular  channel.  Initially,  the  model  was  solved  with  an  assumption  of 
progressive  surface  wave.  It  resulted  in  the  relationship  between  thermocline 
oscillations  on  one  hand  and  amplitude  of  barotropic  tide,  depth  of  the  shelf  and  deep- 
sea  regions,  thermocline  depth  and  density  difference  between  the  bottom  and  surface 
layers  on  the  other.  The  amplitude  ratio  calculated  by  the  model,  using  physical 
parameters  (thermocline  depth  and  density  difference)  corresponding  to  June  2001, 
was  around  10,  which  is  in  agreement  with  the  empirical  results.  However,  the  phase 
difference  between  surface  elevation  and  currents  from  this  model  is  not  appropriate 
for  diurnal  barotropic  tide  in  the  Adriatic,  which  shows  features  similar  to  the 
standing  surface  wave.  Therefore,  a  two-layer  analytical  model  for  the  standing 
surface  wave  was  developed,  in  which  the  Adriatic  was  approximated  by  a  semi- 
enclosed  basin  with  a  step-like  topography.  This  model  covered  two  different  cases. 
In  the  first  classical  boundary  conditions  were  imposed  for  the  surface  and  internal 
mode  at  the  closed  and  open  end  of  the  basin,  giving  standing  wave  both  at  the 
surface  and  at  the  thermocline.  In  the  second  case  radiation  boundary  conditions  were 
imposed  for  the  internal  mode,  resulting  in  progressive  character  of  the  thermocline 
wave.  Both  models  succeeded  in  capturing  the  basics  of  internal  tide  generation,  as  it 
occurs  under  the  influence  of  changes  in  barotropic  flow  at  the  shelf  break.  However, 
they  also  demonstrated  sensitive  dependence  on  the  parameters  that  define  the 
Adriatic  basin  (deep  sea  and  shelf  region  depths,  depth  discontinuity  coordinate, 
position  of  observation  point  in  the  respective  coordinate  system). 
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The  results  from  the  present  report  have  helped  to  define  the  precise  positions  of  the 
ITHACA  measurement  points.  They  are  shown  in  Figure  25,  whereas  the  coordinates 
of  the  ADCP,  thermistor  and  CTD  stations  are  given  in  Table  3.  The  stations  were 
selected  so  as  to  enable  hypotheses  stated  in  this  report  to  be  tested,  but  also  with  the 
aim  of  complementing  stations  that  will  be  occupied  in  the  framework  of  another 
project  planned  for  the  same  area  and  time  (Dynamics  of  the  Adriatic  in  Real  Time  - 
DART).  Let  it  be  added  that  the  instruments  needed  for  the  present  project  (ADCPs, 
thermistors,  CTD  probes,  meteorological  package,  tide  gauge)  are  prepared  for  the 
start  of  the  experimental  work  in  February  2006. 


Figure  25.  Position  of  stations  to  be  occupied  in  the  framework  of  ITHACA  project. 


53 


Table  3.  Basic  information  on  the  ADCP,  thermistor  and  CTD  stations  to  be  occupied 
in  the  framework  of  ITHACA  project. 


Station 

Instruments 

Latitude 

Longitude 

Approx, 
depth  (m) 

ITHACA- 1 

ADCP 

42.64553 

17.14057 

106 

ITHACA-2 

ADCP 

42.74562 

16.68553 

102 

BISEVO 

THERMISTORS 

42.95337 

16.01276 

/ 

SUSAC 

THERMISTORS 

42.75633 

16.51851 

/ 

LASTOVO 

THERMISTORS 

42.72775 

16.90049 

/ 

CTl 

CTD 

43.03689 

16.27861 

95 

CT2 

CTD 

43.00000 

16.33333 

105 

CT3 

CTD 

42.96088 

16.38643 

1 10 

CT4 

CTD 

42.92200 

16.43593 

115 

CT5 

CTD 

42.88311 

16.49072 

115 

CT6 

CTD 

42.84422 

16.54729 

110 

CT7 

CTD 

42.79650 

16.61799 

115 

CT8 

CTD 

42.74562 

16.68553 

102 

CT9 

CTD 

42.69574 

16.75233 

135 

CT10 

CTD 

42.65155 

16.81596 

145 

CTl  1 

CTD 

42.60559 

16.88136 

160 

CTl  2 

CTD 

42.55963 

16.94323 

175 

CTl  3 

CTD 

42.50307 

17.02277 

190 
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